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ABSTRACT 

Context. Our knowledge of the inner structure of embedded massive young stellar objects is still quite limited. Thus, it is difficult to 
decide to what extent the mass accumulation onto forming massive stars differs from the process of low-mass star formation. 
Aims. We attempt to overcome the spatial resolution limitations of conventional thermal infrared imaging. 

Methods. We employed mid-infrared interferometry using the MIDI instrument on the ESO/VLTI facility to investigate M8E-IR, a 
well-known massive young stellar object suspected of containing a circumstellar disk. Spectrally dispersed visibilities in the 8-13 jim 
range have been obtained at seven interferometric baselines. 

Results. We resolve the mid-infrared emission of M8E-IR and find typical sizes of the emission regions of the order of 30 milliarc- 
seconds (= 45 AU). Radiative transfer simulations were performed to interpret the data. The fitting of the spectral energy distribution, 
in combination with the measured visibilities, does not provide evidence for an extended circumstellar disk with sizes ^100 AU but 
requires the presence of an extended envelope. The data are not able to constrain the presence of a small-scale disk in addition to an 
envelope. In either case, the interferometry measurements indicate the existence of a strongly bloated, relatively cool central object, 
possibly tracing the recent accretion history of M8E-IR. In addition, we present 24.5 yum images that clearly distinguish between 
M8E-IR and the neighbouring ultracompact Hn region and which show the cometary-shaped infrared morphology of the latter source. 
Conclusions. Our results show that IR interferometry, combined with radiative transfer modelling, can be a viable tool to reveal 
crucial structure information on embedded massive young stellar objects and to resolve ambiguities arising from fitting the SED. 

Key words, stars: formation - techniques: interferometry, radiative transfer - individual object: M8E-IR 



1. Introduction 

High-mass stars predominantly form in clustered environments 
much farther away from the Sun, on average, than typical well- 
investigated low-mass star-forming regions. Thus, high spatial 
resolution is a prerequisite for making progress in the observa- 
tional study of high-mass star formation. Furthermore, all phases 
prior to the main sequence are obscured by dense circumstellar 
environments. This forces observers of deeply embedded mas- 
sive young stellar objects (MYSOs) to move to the mid-infrared 
(MIR), where the resolution of conventional imaging is limited 
to > 0725 even with 8-m class telescopes. Hence, one traces 
linear scales still several hundred AU in size even for the near- 
est MYSOs, and conclusions on the geometry of the innermost 
circumstellar material remain ambiguous. MIR emission moder- 
ately resolved with single-dish telescopes could even arise from 
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the inner outflow cones (e.g.,|De_Buizer 2006tlLinz et"aLll2005l) . 
A way to overcome the diffraction limit of single telescopes is to 
employ interferometric techniques. We are presently conducting 
a larger survey toward MYSOs based on MIR interferometry. 
While the results for the other objects will be reported in subse- 
quent publications, we concentrat e here on the o bject M8E-IR, 
a prominent BN -type MYSO (cf. )Hennindll990l) at a distance 
of 1.25-1.5 kpc dArias et alJl2007HSimon et al.lll984l) . Although 
M8E-IR had been well investigated in the 1980s, the spatial res- 
olution for most of the I R observations of M 8E-IR was poor. An 
exception is the work bv lSimon et al.l ((1985) who speculated on 
the existence of a small circumstellar disk around M8E-IR based 
on thermal infrared lunar occultation data. Here, we present our 
work on M8E-IR which includes 8-13 fim interferometry to dis- 
sect the MYSO M8E-IR itself, and N- and Q-band imaging to 
have a fresh look at M8E-IR in its relation to the environment. 
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Fig. 1. Left: Q-band image of the M8E region at 24.5 /-im with COMICS/Suba ru. The two pr eviously known objects are annotated. 
The black plus indicates t he position of a 44 GH z methanol maser according to lVarttsI d 1999b . The white plus marks the position of 
the cm continuum source (Molin ari et alj fl998). Right: N-band image of the M8E region at 10.5 /im with COMICS/Subaru. The 
two plus signs mark the same objects as in the left image. The IR counterpart of the cm source is not detected, while the third IR 
source is faintly visible. 



2. Observations 

2.1. Mid-infrared imaging with COMICS at Subaru 

The SMOKA data archive dBaba et all |2002|) contains 8.2-m 
Subaru observations of M8E-IR in the thermal in frared, obtained 
on Ju ne 08, 2004 with the camera COMICS (lOkamoto et al.l 
I2003I) . utilising the 24.5 fim Q-band filter (A c = 24.56 /mi, AA 
= 0.75 jum) as well as the 10.5 yum N-band filter (A c = 10.48 /mi, 
AA = 1.05 /an). For imaging, COMICS employs one 320 x 240 
pixel 2 Si:As IBC detector from Raytheon with a nominal pixel 
scale of 0.13 arcseconds. For the observations, the whole de- 
tector was read out with an elementary exposure time of 0.06 
s. These exposures were repeated 16 times per chop position. 
Chopping with a frequency of 1 Hz, a nominal throw of 20" 
and a chop angle of 140° was applied to remove the strong ther- 
mal background. Since for COMICS/Subaru, the residual pattern 
after chopping subt raction is negligible for most purposes (cf. 
Fuii wara et al] 20061) . additional nodding was not applied during 
the observations. The total integration time (chop-on plus chop- 
off) is 400 s in Q-band and 200 s in N band, respectively. The 
stars HD 148478 (Q-band) and HD 169916 (N-band) were ob- 
served directly before and after M8E-IR for flux calibration. 



2.2. MIDI interferometry at the VLTI 

Visibilities in the mid-infrared wavelength range 8- 13 /mi have 
been obtained with the instrument MIDI dLeinert et al] [2003) 
at the Very Large Telescope Interferometer. Within the frame- 
work of Guaranteed Time Observations for MIDI as well as 
Director's Dicretionary Time, we observed M8E-IR at seven 
baseline length / baseline orientation combinations between June 
2004 and June 2005. In Table Q] we list the UT dates and UT 
times for the fringe track data, the projected baseline lengths 
and the position angles of the projected baselines on the sky 
(counted from north via east on the sky), as well as the tele- 
scope configurations used and the observing proposal numbers. 



Table 1. Log of MIDI observations of M8E-IR. 



UT date and time 


B 


PA 


Telescope 


ESO 






[m] 


[deg] 


pair 


Program 


2004-06-05 


08:07 


96.8 


+42.7 


U1-U3 


073.C-0175(A) 


2004-06-05 


09:57 


82.2 


+44.5 


U1-U3 


073.C-0175(A) 


2004-08-01 


01:51 


46.6 


+38.4 


U2-U3 


273.C-5044(A) 


2005-03-02 


08:58 


46.8 


-85.9 


U3-U4 


074.C-0389(B) 


2005-06-24 


07:37 


51.6 


-42.2 


U3-U4 


075.C-0755(B) 


2005-06-24 


09:29 


43.4 


-10.6 


U3-U4 


075.C-0755(B) 


2005-06-26 


00:26 


55.7 


-06.6 


U1-U2 


075.C-0755(A) 



We refer to Leinert et all (12004 for a more detailed description 
of the standard observing procedure for MIDI observations. For 
all our observations, the so-called HighSens mode was used: dur- 
ing self-fringe tracking, all the incoming thermal infrared signal 
is used for beam combination and fringe tracking, while the pho- 
tometry is subsequently obtained in separate observations. We 
use the MIDI prism as the dispersing element, hence, we finally 
obtain spectrally dispersed visibilities with a spectral resolution 
of R « 30. HD 169916 was used as the main interferometric 
and photometric standard star and always was observed immedi- 
ately after M8E-IR. In addition, all calibrator measurements of a 
night were collected to create an average interferometric transfer 
function and to assign error margins to the measured visibilities. 
For the August 01, 2004 observations, the conditions were al- 
most photometric, the airmass of both M8R-IR and HD 169916 
was minimal (1.01), and the data are hardly affected by atmo- 
spheric disturbances (e.g., ozone feature). Therefore, we use the 
dispersed photometry from this measurement to provide the N- 
band spectrum later used in the SED fitting (see Sect. 13. 3I >. 
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Fig. 2. Left: Collection of dispersed MIDI visibilities. The inset schematically shows the baseline lengths and orientations (cf. Table 
[TJ. Note that the apparently rising slope of the purple UT3-UT4 visibility curve from 24-06-2005 (dashed line) is probably an artefact 
of the measurement at very high airmass (see Sec. 13.2b . Right: The two visibility curves with very similar projected baseline lengths 
(46.6 m vs 46.8 m) but strongly differing position angles (+38.4°, upper curve vs. -85.9°, lower curve). The formal error bars 
(+/— 3(r) are indicated. 



3. Results 

3.1. MIR imaging 

In Fig. [T] we show the 24.5 fim Subaru/COMICS image on the 
left. M8E-IR is still the dominating source at this wavelength. At 
a nominal resolution of 0775, the emission remains compact and 
barely differs from the profile of the supposedly un resolved cal- 
ibrator star. In a very recent MIR survey paper by Ide Wit et al.l 
d2009l) . the same Subaru data set is used to analyse the 24.5 fim 
intensity profile of M8E-IR in more detail and to model the pro- 
file and the SED by means of spherically symmetric model con- 
figurations. These authors come to the same conclusion regard- 
ing the compactness of the M8E-IR profile. A second faint point 
source not yet reported in the literature is visible roughly 6" west 
of it. Furthermore, we clearly detect MIR emission arising from 
the ne i ghbouring radio so urce (Si mon et al.ll984tlMolinari et alj 
ll998h . lSimon etalj d!985l) had already reported a detection of 
this source in N- and Q-band (fluxes only, measured with a 6" 
diaphragm). For the first time, the Subaru MIR imaging spatially 
resolves this emission. It is comet-shaped, with the apex directed 
away from M8E-IR. This morphology could be an intrinsic prop- 
erty of this UCHu region, or it could be sh aped by the molecu- 
lar ou tflow probably arising f rom M8E-IR dMitchell et alJI 19881 
1 19921) onto this radio source. Damian i et al.l d2006l) reported an 
X-ray source * 1'.'7 southwest of M8E-IR ([DFM2004] 845). 
Accounting for the uncertainty of the X-ray position (2'.'56), the 
X-ray emission can still be associated with M8E-IR, but not with 
the known radio source or the abovementioned faint MIR point 
source. 

At 10.5 fim we easily detect M8E-IR (Fig. [T] right). The source 
is very bright, which causes some additional image artefacts 
(diffraction spikes in diagonal directions as well as a detec- 
tor "drooping" effect, causing multiple fainter copies of the 
strong source along the north-south direction). Hence, we will 
not further investigate the shape of M8E-IR in this image. 
The radio source counterpart is not detected at 10.5 fim (rms 
noise * 8 mJy). Considering the COMICS filter characteristics 
(S ect. 12. IK combined w ith the detection of the radio counterpart 
by ISimon et all d 1985b in a broad N-band filter (0.97 Jy), this 
implies a very strong and broad 9.7 fim silicate absorption fea- 



ture or time variability. The third source in the 24.5 fim image is 
detected at a 4cr level at 10.5 fim. 

3.2. MIR interferometry 

We reduced the interferometric data with the MIA+EWS pack- 
age, version 1.5, developed at the MPIA Heidelberg and the 
University of Leiden. The resulting visibility curves are col- 
lected in Fig. |2] left. The object is clearly resolved in all con- 
figurations with visibilities between 0.09-0.35. If we assume a 
Gaussian intensity distributiorQ of the source, the visibilities in- 
dicate an intensity FWHM of « 20-25 mas (8.5 fim) and 32- 
38 mas (12.0 fim) which is in rou gh agreement with the exten- 
sion of the small component of Sim on et all (fl98 5). We note 
that these visibilities, although not rea ching the relatively high 
levels of most Herbig Ae/Be stars (e.g. jLeinert e t al. 2004), are 
qualitatively different from the very low visibilities (0.01-0.05) 
found for several of the other objects in our samp le as well 
as recently reported for two other massive YSOs dde Wit et al.l 
120071: IVehoffet al.1l2008l) . Still, to learn more about the poten- 
tially more complicated intensity structure, a simple Gaussian 
is not a sufficient ansatz, especially if the object is strongly re- 
solved by the interferometer. Further modelling is necessary for 
interpretation. 

Can we infer the geometry of the source from our visibil- 
ity curves alone? Most visibilities have been taken at differ- 
ent baseline length / position angle combinations. Still, there 
are two measurements with almost the same projected base- 
line length (46.6 m vs 46.8 m), but position angles differing by 
55.7° (cf. Table [T). A spherically symmetric intensity distribu- 
tion will result in visibilities not being a function of position 
angle. These two measurements, however, show different visi- 
bility levels (Fig. [2] right). But when considering the 3-cr error 
bars, this difference is not fully conclusive, since the error bars, 
at least in the 8 to 1 1 fim range, still overlap slightly. At longer 
wavelengths, the two error intervals are almost distinct. While 



1 The relation between the visibility and the intensity full width half 
maximum is given by V(u) = exp(-7r 2 /(4 ln2) u 2 FWHM 2 ) in this case, 
where the spatial frequency is u = A/B, with B being the modulus of 
the projected baseline length on the sky. 
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Table 2. Parameters of the best-fitting Robitaille model 
3007526. 
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Fig. 3. SED of M8E-IR. Crosses = observed fluxes, diamonds = 
upper limits. The solid line gives the best-fitting radiative trans- 
fer mo del 3007526 from the grid of models by Robitaille et al. 
(2007) (effective synthetic aperture ~ 2"). The dashed black 
curve shows the SED of the bloated central star used in that 
model. The red dotted line shows the spherical model with a hot 
star (also mentioned in Sec. [33} that fits the SED but not the visi- 
bilities (cf. Fig. [5]l. The inset shows a zoom into the 10/mi region 
using a synthetic aperture < 1", appropriate for the underlying 
MIDI spectroscopy data (shown as plus signs). For legibility, the 
symbol sizes are similar to or larger than the real error bars. 



we cannot draw strong conclusions from these findings, this is a 
first hint that the mid-infrared intensity distribution of M8E-IR 
is indeed not fully spherically symmetric. 
Visibilities can be considered as the ratio of correlated flux to to- 
tal flux. The N band total flux spectrum of M8E-IR shows a mod- 
erately deep silicate absorption spectrum (see inset of Fig. |3j. 
Most of the curves displayed in Fig. [2] (left) show a drop in the 
visibility level over the silicate feature. However, the visibility 
curve corresponding to the 2005-06-24 09:29 UT measurement 
does not show this behaviour, but instead rises monotonically 
from 8 to 13 /mi. We think that this is an artefact of the observa- 
tion. While all the other measurements were taken at intermedi- 
ate elevations (airmass < 1.73), the data set in question was ob- 
tained at a very low elevation (airmass 2.41) in order to get this 
(u,v) point, corresponding to a relatively short projected base- 
line length. At such a low elevation, the performance of MIDI 
is certainly affected. Since the adaptive optics system cannot be 
locked on M8E-IR itself, but only on a faint optical star more 
than 40" away, the beam stability and therefore the beam over- 
lap on the MIDI detector may be compromised. Since the PSF 
is smaller at 8 micron than at 13 micron, a smaller beam overlap 
will lead to apparent losses in correlated flux preferentially on 
the short-wave side of the spectrum. We have critically checked 
the overlap for this measurement and indeed find that the two 
beams are more separated than for other observations of M8E- 
IR as well as for the calibrators. By artificially restricting the 
overlap mask to a smaller central area, the slope of the visibility 
curve can be flattened with (calibrated) visibility values at 8 mi- 
cron slightly higher than at 13 micron. Hence, the low 8 micron 
correlated flux of the indicated visibility curve is an artefact of 
the measurement process. Therefore, we will not further inter- 
pret the shape of this particular visibility curve. 



Parameter 


Value 


Parameter 


Value 


Effect, temperature 
Stellar mass 
Outer disk radius 
Scale height exp.' 1 ' 
Envelope radius 
Ay along LOS 
Cavity cone angle 


4,740 K 
13.5 M Q 
16 AU 
1.05 
6.8E+4 AU 
33 mag 
3.3° 


Stellar radius 
Disk mass 
Density exp.' 1 ' 
Viscosity a par. (1) 
Mass infall rate <2) 
Inclination angle 


125 Re 
0.71 M Q 
2.05 
1.3E-2 
7E-5 M Q /yr 
18° 


(1) for the disk <2) 


"or the envelope 


Dust model: Ki 


m et all d 1994b 



3.3. Modelling 

We apply self-consistent continuum radiative transfer modelling 
to M8E-IR in order to produce synthetic MIR intensity maps and 
to compare their spatial frequency spectrum with the observed 
visibilities. Self-consistency here means that the dust tempera- 
ture is calculated from the condition of (radiative) energy conser- 
vation. Hence, the luminosity sources in the system (central ob- 
ject + accretion) initially heat up the dust until a radiative equi- 
librium is reached, the temperature (for the given dust density 
structure) conve rges, and the total luminosity is recovered in the 
emerging SED (iBiorkman & Woodll200U IWhitnev et al.ll2003l: 
Men'shch ikov & Henningll997l) . We are mainly concerned with 
the question of which spatial distribution of the circumstellar 
material can account for both the SED and the visibilities of 
M8E-IR. 

F or SED fitting, we us e the M8E-IR photometric data collected 
inlMueller et alj d2002l) plus new 1 .2 mm data from Beltran et al. 
d2006l) . We stress that no (sub-)millimeter interferometry on 
M8E-IR is reported in the literature which could spatially disen- 
tangle the flux contributions from M8E-IR and the radio source 
8" away. Hence, we consider the fluxes for A > 24.5 fim as up- 
per limits in our modelling. Furthermore, we in clude new optical 
photom etry in the B, V, and I filters reported by Prisi nzano et al.l 
( 120051 their object CI* NGC 6530 WFI 13458) as well as our 
new 24.5 /mi photometry (220 + 20 Jy). The 8-13 /mi total flux 
spectrum taken during the MIDI measurements (Sect. 12.21 ) was 
used to further curtail the number of via ble models. 
We used the SED online fitting tool of iRobitaille et alj d2007l) 
that can in principle utilise models including an envelope plus 
circumstellar disk. We refer to this publication for details on the 
setup of these models. To produce the synthetic MIR intensity 
maps corresponding to the Robitaille models, we use the un- 
derlyin g Monte Carlo radiative transfer code of I Whitney et al.l 
d2003l) . The SED fitting tool delivers the full parameter set 
for the best-fitting models which can be directly included 
in the Whitney code. The resulting images are then Fourier- 
transformed, and cuts through this spatial frequency spectrum 
of the images give the synthetic visibilities. 
Intriguingly, we do not find a unique solution from this SED fit- 
ting. The formally best SED fit points to a model featuring a 
very compact circumstellar disk (< 20 AU), a larger surround- 
ing envelope with small bipolar cavities and a bloated cooler 
central object (cf. Table |2]i In Fig. [3] we show this best SED fit 
to M8E-IR by the Robitaille models as black continuous line. 
Due to the linkage of the Robitaille grid to evolutionary tracks, 
certai n size parameter s cannot be independently chosen by the 
fitter (Robitaille 2008). This especially affects the possible disk 
size which is not well constrained in our case. The data would 
still allow for a somewhat larger disk. We have repeated the cal- 
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culation of the intensity maps with all model parameters fixed to 
the values of the model best fitting the SED, but have used larger 
disk sizes of 50, 100, 250, and 500 AU. For disk radii ^ 100 AU, 
the resulting MIR visibilities eventually drop. Hence, the differ- 
ence to the values measured with MIDI increases, which makes 
such solutions less likely. 

Among the well-fitting Robitaille SED models are also configu- 
rations without a compact disk (axisymmetric flattened envelope 
+ outflow cavities only), which nevertheless give almost identi- 
cal elevated visibilities (see Fig. [5]). This suggests that in the case 
of M8E IR, the choice of the central object actually governs the 
resulting visibility levels (see below). 

We have also compared the spectral shape of the visibility curves 
with the models. For the comparison, we have computed syn- 
thetic images based on the parameters from the SED fits at 
eleven wavelengths between 8 and 13 micron. Their Fourier 
transform is then probed at spatial frequencies according to the 
projected baseline lengths of the MIDI measurements. In gen- 
eral, the Robitaille models predict somewhat lower visibilities 
than indicated by the measurements. The general shape of all 
the synthetic visibility curves is similar. Therefore, we show as 
one example the 46.8-m baseline comparison in Fig. [4] We com- 
pare the best envelope + disk configuration (Robitaille model 
3007526, see Table |2]i and the best envelope-only configuration 
(Robitaille model 3004120), both having a cool central object 
(see also Fig- [5j, to the MIDI measurement. The model visibil- 
ity curves show a clear decrease in the visibility level toward 
the amorphous silicate feature. The modelled shape would arise 
naturally in the circumstellar disk context when the correlated 
flux comes mainly from an inner, optically thicker region with 
a somewhat suppressed silicate emission feature, while the total 
flux is domi nated by more optically thin emission from the wider 
disk surface (van Boekel et al .120051) . However, as seen in Fig.|4] 
the envelope-only model also gives the same behaviour. Hence, 
based on the modelled visibilities, it is not possible to distin- 
guish between the two cases for our special setup. This is under- 
standable considering the compactness of the disk in question, 
which would not be resolved even by our longest baseline mea- 
surement. The envelope by far dominates the mid-infrared emis- 
sion of the whole object. The measured visibilities also show a 
level decrease. However, it is not centered around 9.7 micron, 
but occurs at longer wavelengths. Currently, we do not under- 
stand this shift. It is not clear if it can be fully explained by 
differing dust properties in M8E-IR and assumed for the radia- 
tive transfer simulations (opacities an d grain size model from 
lLaor&Drainelll993t iKim et alj[l994b . To further analyse and 
reproduce the spectral shape details of the visibility curves is be- 
yond the scope of this paper. 

We explicitly mention the low inclination angle of < 30° in- 
herent in all our Robitaille models fit ting the SEP data. This 
is in contrast to the disk hypothesis of ISimon et alj (1 19851) that 
featured a very large inclination. However, Simon et al. saw a 
more or less symmetric intensity distribution in their 10 fim lu- 
nar occultation data; hints for deviations came mainly from lunar 
occultations recorded at 3.8 //m. At this shorter wavelength, con- 
figurations with disk + envelope including outflow cones are nat- 
urally more structured, and scattered light in the inner regions of 
these cones can still contribute, in particular if grains larger than 
the typical 0. 1 /im dust particles are involved. 



Comparison MIDI vs models (46.8 m baseline) 
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Fig. 4. Spectral shape comparison of the visibilities for the 46.8- 
m baseline as an example. The MIDI measurement plus 3-cr er- 
ror bars is given as a solid line. Furthermore, radiative trans- 
fer synthetic images between 8 and 13 micron have been pro- 
duced for the best envelope + disk model (dashed line) and the 
best envelope-only model (dotted line). Cuts through their spa- 
tial frequency spectrum corresponding to the 46.8-mbaseline are 
concatenated to display the synthetic spectral visibility shape of 
these models. 

3.3.1. Comparison with canonical 1D models 

As mentioned above, the Robitaille models are not optimally 
sampled in the parameter ranges appropriate for massive YSOs. 
Formally, the Robitaille models cover a sufficiently large param- 
eter space regarding the 14 important parameters. In particular, 
hot, non-bloated central stars are includ ed, up to effective tem - 
peratures of > 50,000 K (see Fig. 3 in IRobitaille et al.ll2006l) . 
However, not all parameter combinations are available for the 
fitting since it is not possible to explore the huge parameter 
space in its entirety in a completely unbiased manner when con- 
structing the model grid. Hence, the set of models available for 
the SED fitting does not form a complete "grid" in a technical 
sense, since paramet ers are randomly sampled within ranges (cf. 
IRobitaille et all 120061 especially the sections 2.1.4 and 2.2.1). 
Since we have found quite peculiar models from our fitting 
attempt (see above), we felt compelled to do a cross-check of 
our results. We could have missed simpler models among the 
top-rank fitting results that nevertheless agree with the observed 
SED. Traditionally, such BN-type objects have been modelled as 
a sph erical dust shell surrounding a hot, non -bloated central star 
(e.g., iGuertler et alJl!991t Ide Wit et alJl2009h . To check if even 
such canonical configurations can account for the SED and the 
visibiliti es, we employed well-tested mod el s, based on the code 
used b y [Men'shchikov & Henningl l[l997); Men 'shchikov et all 
(1999). Unlike the Whitney code, the Men'shchikov program 
is a grid-based radiative transfer code. It comes with its own 
ray-tracing module to produce high-definition intensity maps. 
Here, we can tune the parameters without being restricted by 
a predefined parameter grid, and the fine-tuning process in the 
case of pure ID models can be much faster. With a spherically 
symmetric envelope (density power law ~ r~ 1J , M env = 40 M ) 
and a 30,000 K central star with 3.7 R H , it is possible to fit the 
SED of M8E-IR well (Fig. [3] see also Ide Wit et all d2009lPI) . 



2 These authors find slightly different model parameters, because 
they try to find a balance between fitting the 24.5 /im intensity pro- 
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Still, compared to the measured visibilities, these models result 
in far too low visibilities (< 0.05) for M8E-IR in the baseline 
range 30-60 m over the whole 8-13 /mi range, even when we 
add a population of very large grains (200 yum) to the central 
shells to raise the visibilities. The error margins of the MIDI 
visibilities (on average 10%) do not account for such large 
differences. Independently, we found that the inclusion of a cold 
(~ 2, 000 K) central object will sufficiently raise the visibilities 
also in the Men'shchikov models. The circumstances are illus- 
trated in Fig. [5] where the (u,v)-spectra of the 12 fim synthetic 
images based on 3 different models are compared. At such a 
longer wavelength outside the silicate absorption feature, the 
influence of different dust models on the results should be small. 
Furthermore, potential observational biases due to decreased 
beam overlap during the MIDI measurements are smaller at 
longer wavelengths (see Sec. |3.2| >. Although the (u,v)-spectra of 
the cool star models still show somewhat lower visibilities than 
the measured ones (the visibilities are not included in the fitting 
loop), obviously they are qualitatively different from the hot star 
model. 



3.3.2. Feasibility of a cool central object 

Interestingly, the best-fitting models for M8E-IR in the Ro- 
bitaille model grid show central stars of 10-15 M that are 
strongly bloated and, therefore, have relatively low effective 
temperatures (< 5000 K). Such solutions can occur since the 
Robitaille grid also comprises the full range of canonical pre- 
main sequence evolutionary tracks from the Geneva grourJ3 
as possible parameterisation of the central objects. Although 
M8E-IR cannot straightforwardly be identified with correspond- 
ing very early evolutionary stages, the tendency for a bloated 
central sta r in the case of M8E-IR may be valid. As found 
already by Kippenhahn & M ever-Hofmeisterl dl977l) . accretion 
with high rates onto main sequence stars can temporarily puff 
up su ch stars. Furthermore, iHosokawa & Omukail (120081 [2009) 
and lYorke & Bodenheimerl (120081) have recently computed the 
pre-main-sequence evolution of stars as a function of the accre- 
tion rate onto the forming star. They find that for accretion rates 
reaching 10~ 3 M o /yr, the protostellar radius can temporarily in- 
cre ase to > 100 R B , in ac cordance with our model fitting. Note 
that lMitchell et al.f(ll988l) revealed high- velocity molecular out- 
flows from M8E-IR based on M-band CO absorption spectra and 
speculated on recent (< 120 yr) FU Ori-type outbursts for this 
object. If these multiple outflow components really trace recent 
strong accretion events, the central star indeed could have been 
affected. Furthermore, M8E-I R is not detected by cm observa- 
tions with medium sensitivity (ISimon et al.ll984_lMoIinari et al.l 
1998). This could be explained by large accretion rates still 
quenc hing a forming hypercompact Hn region (e.g., Walmsley 
1 19951) . Similarly, a bloated central star with T eS <c 10000 K also 
gives a natural explanation of these findings. 



file and properly fitting the SED. Thereby, they take some of the exist- 
ing single-dish (sub-)millimeter points as face value. There are indica- 
tions, however, that at long wavelengths the neighbouring UCHu region 
(Fig.0 left) dominates the emission of the region. Therefore, these long 
wavelengths fluxes should be considered as strict upper limits. 
3 http://obswww.unige.ch/~mowlavi/evol/stev_database.html 



RT models vs observations for M8E-IR at 12 micron 




Baseline [meter] 

Fig. 5. Cuts through the 12 fj.m spatial frequency spectrum (in 
units of the interferometric baseline length) of 3 models well fit- 
ting the SED. The black solid line refers to a model with a small 
circumstellar disk, flattened envelope and cool bloated central 
object (Table [2}. The blue dashed line denotes a model with a 
flattened envelope and cool bloated central object without a disk; 
the intensity distribution is almost identical to the first model. 
Due to the low inclination of these 2D models (18°), the cen- 
tral intensity distribution shows only small deviations from ra- 
dial symmetry. Hence, differences in the model visibilities with 
varying position angle remain small (< 10%), and the cuts can 
be adopted as representative. The classical approach (hot central 
star + spherical envelope, red dotted line), however, gives far 
too low visibilities compared to the MIDI data (plotted as plus 
signes including the formal 3<x error bars). 

4. Conclusions 

We have observed the BN-type massive young stellar object 
M8E-IR with the MIDI MIR interferometer at the VLTI. We 
find substructures with MIR sizes around 30 mas. The mea- 
sured elevated visibilities indicate that the usual approach for 
BN-type objects (a spherically symmetric envelope with a hot 
central star) fails in the case of M8E-IR. Most interestingly, our 
data are consistent with M8E-IR harbouring a 10-15 M central 
star that has been bloated (R > 100 R Q ) by recent strong accre- 
tion events. However, with the present data we cannot clearly 
infer the existence of a circumstellar disk since disk-envelope 
and envelope-only models result in very similar SEDs and visi- 
bilities. If a disk exists it cannot be large (;$> 100 AU). Until (sub- 
millimeter observations with sufficient spatial resolution con- 
strain the long-wavelength emission, disks with radii from 15- 
100 AU seem probable. Both aspects, the central star bloating 
and the existence of a circumstellar disk, should be the topic of 
future investigations. For instance, by utilising high-resolution 
IR spectroscopy one can test whether hydrogen recombination 
lines show imprints of di sk rotation along the guidelines of the 
disk-wind simulations bv lSim etaL (2005b. 
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